Non-Entangling Channels for Multiple Collisions of Quantum Wave Packets 
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We consider multiple collisions of quantum wave packets in one dimension. The system under 
investigation consists of an impenetrable wall and of two hard-core particles with very different 
masses. The lighter particle bounces between the heavier one and the wall. Both particles are ini- 
tially represented by narrow Gaussian wave packets. A complete analytical solution of this problem 
is presented on the basis of a new method. The idea of the method is to decompose the two-particle 
wave function into a continuous superposition of terms (channels), such that the multiple colli- 
sions within each channel do not lead to subsequent entanglement between the two particles. For 
each channel, the time evolution of the two-particle wave function is completely determined by the 
motion of the corresponding classical point-like particles; therefore the whole quantum problem is 
reduced to a classical calculation. The calculation based on the above method reveals the following 
unexpected result: The entanglement between the two particles first increases with time due to the 
collisions, but then it begins to decrease, disappearing completely when the light particle becomes 
too slow to catch up with the heavy one. 

PACS numbers: 03.65.Nk,03.65.Yz,03.67.Bg 



In quantum mechanics, particles are often represented 
by wave packets. According to the Schrdinger equation, 
the wave packet of a free particle typically spreads with- 
out a limit. This spread was of concern in the early 
days of quantum mechanics for Schrdinger himself [lj, 
because it had never been observed for macroscopic sys- 
tems. He was able to identify the non-spreading coherent 
wave packets only for the harmonic potential. In the free 
particle case, the spread of the wave packets and the ac- 
companying winding of the quantum phase significantly 
complicate practical calculations, in particular, in the 
treatment of multiple collisions between particles. Typ- 
ically, each collision between two quantum wave packets 
leads to the entanglement between them HH. As a re- 
sult, the wave function of the system after many collisions 
becomes increasingly intractable. 

The studies of colliding wave packets are often moti- 
vated either by the agenda of the controlled generation of 
entanglement Oa, [5|-|8| or by the studies of quantum deco- 
herence 0, 0, Q . In the both contexts, the treatment of 
multiple collisions between quantum wave packets is of 
significant interest. In the former context, it was investi- 
gated for two cold atoms in a harmonic trap in Ref.Q. 

In this article, we consider the problem of a light parti- 
cle bouncing between a heavy particle and a wall, shown 
in Fig. [I] and propose a new method for dealing with mul- 
tiple scattering of quantum wave packets. This method is 
based on the observation that there exist special initial 
conditions which do not lead to entanglement between 
two colliding Gaussian wave packets. We refer to states 
with these conditions as "non-entangling channels" . The 
method consists of decomposing the initial wave function 
of the system as a superposition of non-entangling chan- 
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FIG. 1. (Color online) The system considered: the light par- 
ticle bouncing between the heavy particle and the wall. 



nels. The calculation for each channel becomes largely 
classical, even though the complete solution retains all 
the usual quantum features. For the two-particle prob- 
lem considered, we report a remarkable result: Even 
though the wave packets of the two particles become en- 
tangled after initial collisions, the entanglement begins 
to decrease later and eventually disappears completely. 

The setting of our problem is shown in Fig. [TJ We 
use the variable x for the coordinate and for all sub- 
scripts referring to the light particle. Correspondingly, 
the variable y refers to the heavy particle. The two par- 
ticles have very different masses m x and m y such that 
e = ^Jm x /m y <C 1. We assume the particles to be non- 
relativistic. They interact via the hard-core potential. 
The wall is represented by an infinitely high potential 
step. The Hamiltonian of the problem is 



H = 



Pi 
2rn T 



Py 
2m, 



+ S 1 Q(x-y) + S 2 e(-x), (1) 
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where p x and p y are momentum operators and 0(x) is 
the Heaviside-step-function. We only consider the "hard- 
core" limits Si — > oo and S% — > oo, therefore the light 
particle always remains localized between the heavy par- 
ticle and the wall. 

Initially, the light particle is moving while the heavy 
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particle is at rest. We assume a nearly factorized initial 
state of two Gaussian wave packets (with h = 1): 

Vo (x, y ) ~ exp < — ^ h ip x0 x — ^ > , (2) 

where < xmo < Vmo are the maxima of the Gaus- 
sians, p X Q > the initial momentum of the light particle 
and <jq x and o§ y the widths of the respective wave pack- 
ets, which are assumed to be much smaller than any of 
the three distances: imo, Umo and j/mo — imo- The sign 
"~" in Eq.® means that the wave function needs to be 
slightly modified to remove the overlap of the exponen- 
tially small tails of the particles' wave functions and the 
penetration of these tails into the wall. A possible form 
this modification is given in We further consider 

only sufficiently large values of p x o, such that the widths 
of both wave packets remain much smaller than j/mo until 
the collisions between the two particles stop (see below) . 
Therefore, we use the approximation that every collision 
involves only two constituents: either the two particles 
or the light particle and the wall. 

Let us first recall that a Gaussian wave packet, initially 
with the maximum at smO) the width oq x and momen- 
tum p X Q evolves according to a free-particle Schrdinger 
equation as 

d>(x; t) = N(t) exp j- + w} , (3) 

where (3 x (t) 2 = On x ( 1 H ^3- J is the width param- 

eter, x c (t) = £mo + is the time-dependent posi- 
tion of the maximum and N(t) is the normalization fac- 
tor. The probability distribution of positions |0(x;t)| 2 
remains Gaussian. Its width grows with time as \/3 x (t) \. 
We also introduce the analogous width parameter for the 

heavy particle f3 v (t) 2 = cr2 (H The above- 

y \ rn y (7 0y J 

mentioned narrow- wave-packet assumption implies that 

IA/(*)l.lflr(t)|«i/MO. (4) 

We use the description of the reflection of the light 
particle from the wall given in Ref.0. This description 
incorporates the wall by considering the antisymmetric 
superposition of the free-particle Gaussian and its mir- 
ror image with respect to the position of the wall, i.e. 
<f>(x;t) — (f)(—x;t). Far away from the wall <p(x;t) domi- 
nates before the collision, while (f>(—x; t) dominates after 
the collision. Hence, the description of the particle-wall 
reflection amounts to the transformation x — > — x in the 
free-particle Gaussian (lS3| fiof . 

The single collision of two wave packets via a hard- 
core potential in one dimension has been described in 
Refs.Jl, 0, @|- If the Hamiltonian of the two-particle in- 
teraction is rewritten using the center-of-mass coordinate 
R = (m x x + m y y) / \m x + m y ) and the relative coordi- 
nate r = x — y, the problem separates into a free-particle 
motion for variable R and the wall reflection problem for 



variable r. By analogy with the particle-wall collision, 
the two-particle collision can be described by the trans- 
formation r — > — r with the subsequent transformation 
back to the variables (x,y) If the initial wave func- 
tion has the form ((2J, then, after one collision between 
the two particles, the wave function becomes 

Mx,V,t) ~ exp(Af x x 2 +A v 1 y y 2 + Al y xy+Bfx + Bfy}, 

where all coefficients are functions of time and other pa- 
rameters of the problem. In general, this wave function 
is entangled in terms of variables x and y because of the 
term A* y xy amsearing in the exponent. However, this 
term vanishes [2 0, i.e. A\ v = 0, when[ll| 

m x a$ x = m y a 2 y . (5) 

In this case, not only the entanglement vanishes but also 
the width parameters j3 x {t) and P y (t) continue evolving 
as in the respective free-particle cases, i.e. after the col- 
lision the wave function is 

, / .\ J [x-X U (t)f . . 1 fa y. 

ipi{x,y;t) ~exp^ — — +ip x ix> (6) 

where p x \ and p y \ are the momenta and Xyi(t), Viiif) 
are the maxima of the wave packets moving as the coor- 
dinates of two reflected classical particles. The relation 
§5§ implies that m x (3 x (t) 2 = m y (3 y (t) 2 holds both before 
and after the collision. It also implies that the kinetic 
energies stored in the momentum distributions of each 
of the two particles in the rest frames of their respec- 
tive wave packets are equal to each other - the situation 
somewhat reminiscent of the thermal equilibrium. 

We observe further: Since the evolutions of the pa- 
rameters p x (t) and j3 y (t) remain unaffected after a two- 
particle collision in the case of the condition ([5]) and the 
same is true for the particle-wall collision, the ansatz of 
the form ((6]) remains valid after multiple collisions of the 
light particle with the wall and with the heavy particle. 
The evolutions of the centers of the Gaussian wave pack- 
ets xm (t) and j/m (t) and their respective momenta can be 
calculated by identifying the wave packet's maxima with 
the coordinates of classical point-like particles. This is 
what we call "non-entangling channels" . 

We now turn to the general case when the wave packets 
do not obey the relation ([S]). We limit ourselves to the 
case 

m x al x < m y al y . (7) 

(The opposite inequality affords qualitatively the same 
treatment.) The inequality (J7]) is equivalent to the state- 
ment that the initial wave packet for the heavy particle is 
too broad to obey the relation ([5]) . Central to our method 
is the representation of the initial wave packet for the 
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heavy particle as a superposition of Gaussian wave pack- 
ets with different maxima but the same width, which is 
smaller than the initial one ao y such that it obeys Eq. ([5]) . 
That is, each of the wave packets in this decomposition 
together with the wave packet for the light particle de- 
scribe a non-entangling channel. In this way, we reduce 
the whole quantum problem to the classical calculation of 
the probability distribution for the non-entangling chan- 
nels. 

We divide a^ y into two parts 



'Oy — A(T y,Q + °yTi ( 8 ) 

where a 2 T is defined by the condition m x af )x = m y <jy T 
[cf. Eq.([5])]. and rewrite the initial wave function in 
Eq.Q using the mathematical fact that the convolution 
of two Gaussian functions is a Gaussian function again: 

Xmo} 2 



i>o(x,y) ~ exp 



2^L 



WxOX 



(9) 



dy m exp^ 



bin 



VM0\ 



2A <o 



exp< 



The first Gaussian function inside the integral can be 
interpreted as the distribution of the maxima y m for the 
heavy particle in different non-entangling channels. We 
also introduce the variable x m for the maxima of the light 
particle for each non-entangling channel. According to 
Eq.©, all x m initially coincide with xmo- 

The problem is now reduced to calculating the time 
evolution of the positions x m and y m of classical particles 
and then obtaining the probability distribution of these 
positions, which we denote as P X y,n(^m, ym', t). Here in- 
dex n represents the number of collisions between the 
two particles that occurred before time t in each non- 
entangling channel. Given our assumption of sufficiently 
narrow wave packets, we only consider the moments of 
time between the collisions when particles in each non- 
entangling channel have experienced the same number of 
collisions and, otherwise, are sufficiently far from each 
other and from the wall. Therefore, assigning the same 
index n to all non-entangling channels at a given mo- 
ment of time is justified. Since the momenta of the 
particles do not change between the collisions, their re- 
spective values p xn and p yn only depend on the num- 
ber of preceding collisions n but not on the channel and 
not on time. Below we also use marginal probability 
distributions P^ n (x m ;t) = dj/ m P xy ,„(x m , y m ; t) and 

The time evolution of the wave function of the system 
can now be represented as 



ipn(x,y;t) 



t) (10) 

— oo J — oo 

, [x~x m ] 2 . 1 / [y-y m ] 2 , . 1 

where sign ± implies +, when the light particle moves 
away from the wall, and — , when it moves towards 
the wall. 



According to Eq. (j9j) , the initial probability distribution 
for the underlying classical problem is 



-Pxy,o(^m, y m ; 0) = N Q exp • 



2/moJ 



2A <o 



(11) 



x<5 (x m - xmo) , 



where Nq is the normalization constant. The correspond- 
ing marginal probability distributions P x ,o(a;m;0) and 
-Py,o(2/m;0) are illustrated in Fig. [2j The light particle 
has a definite initial position xmo and the momentum p x0 , 
while the heavy particle has a Gaussian distribution of 
possible initial positions y m and no momentum. The sub- 
sequent evolution of P X y,n(x m , y m ; t) is obtained below in 
the leading order in the small parameter e = y/m x /m y . 

In each two-particle collision, the light particle trans- 
fers momentum to the heavy particle. Eventually, the 
collisions stop when the light particle becomes too slow 
to catch up with the heavy particle. The last collision 
will have index n max ~ j- e [lQj | . 

The evolution of the probability distribution 

-fxy,n(^-m? 

y m ;t) between the collisions, i.e. for the 
fixed value of n, just amounts to the shift of the values 
of x m and y m associated with the motion with respective 
velocities p xn /n^x and p yn /m y . Therefore, the shape of 
(xm, y m ; t) only changes when n changes. 
Two trends in the evolution of the probability distribu- 
tion Pxy ,n{X-mi 2/m 

; t) as a function of n can be anticipated, 
namely, the initial increase in the spread of the possible 
values of x m and the initial decrease in the spread of 
the possible values of y m . Both trends originate from 
the time delay between collisions for different initial po- 
sitions of the heavy particle y m o'- The larger y m o, the 
later the first and the subsequent collisions occur. The 
effect of each such time delay is rather small, but it ac- 
cumulates after many collisions. The resulting evolution 
of the particle coordinates becomes (Toj |: 

/ \ / \ sin(2en) 

Xm{ymO;t) = X M {t) H (2/mO-?/Mo), (12) 

e 

y m (ymo;0 = 2/m(*) + cos(2en)(y m0 - Vmo), (13) 

where XM(t) and yyi{t) represent the "reference" trajec- 
tory for the center of the probability distribution corre- 
sponding to y m0 = yuo, and x m (y m0 ;t) and y m (ymo;t) 
stand for the positions of the two classical particles at 



m. 



Px,o(z m :0) 








FIG. 2. (Color online) Classical probability distributions: 
green - P x , (x m ; 0); red - P y fi(y m ;0). 
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time t provided that the heavy particle was at y m o at 
t = 0. The above equations reveal that the distances 
Xm(Vmo; t)-x M (t) and y m (y m o;t)-yu(t) are proportional 
to the initial distance y m (y m o; 0)-y M (0) = y m o-yuo and 
the scaling factors do not depend on y m o; in particular, 

x m (y m0 ;t)-x M (t) = tan ( 2e ") [y m (y m0 ;t)-y M (t)}. (14) 



Given the initial distribution P X y, 0(^111; ym', 0) and the 
scaling relations (flllf . (IT4")) . we obtain the distribution at 
later times as 



^xy,n(^mi ym, ^) — N n exp 



[ijm - 2/M(t)] 



2Aa?, 



(15) 



y,n 



. . . . tan(2e?i) , 
xd £ m - XM(t) [y m - j/m(*)J , 



where iV„ are normalization factors and 
^&y,n = Acy.ol cos(2en)| 



(16) 



[cf. Eq. (|13p ]. The corresponding marginal distribution 
Py,n(yml t) is, up to a prefactor, equal to the exponential 
factor in Eq. (|15|) . while the marginal distribution with 
respect to x m is 



P x ,n(x m ) = exp ^ - 



[x m - x M (t)} 



(17) 



where Aa x<n = A<J y ,o\ sin(2en)|/e [cf. Eq. (fT2)) ]. 

According to Eqs. ([l~5"T) and (p~6|) . the marginal distribu- 
tion Py, n (ym',t) contracts to a (5-function when the num- 
ber of collisions reaches the critical value n cr « ^ . Since 
this value is equal to n max , the collisions stop at the same 
time. 

According to Eq. flTUl) . there is a direct connection be- 
tween the entanglement properties of the wave func- 
tion of the system in terms of variables x and y and 
the factorizability of the classical probability distribution 
y m ;t) in terms of variables x m and y m . The 
former is entangled, if and only if the latter is not factor- 
izable. For example, the initial wave function ipo(x,y) is 
not entangled and the corresponding probability distri- 
bution P X y, 0(^111, ym', 0) given by Eq. lfTTj) is factorizable. 
According to Eq. ([T5|) . the distribution P xy , n (x m , y m ; t) 



becomes not factorizable at later times due to the pres- 
ence of y m in the 5-function. Therefore, the wave function 
of the system becomes entangled. However, after the crit- 
ical number of collisions n cr , the Gaussian in Eq. (|15l) con- 
tracts to the (5-function and as a result P Ky , n (xni, ym', t) 
becomes factorizable again. Thus, the two particles are 
no longer entangled at the time when the collisions stop. 

The narrow- wave-packet assumption Q remains valid 
until the last collision if (em x ao x v x o) / (jrh) 3> lflCj- 

Substituting Eq.flTSl) into Eq. (fTU|) . we obtain 



Mx,y,t) ~ exp{A x n x x 2 +Aiy y 2 +A*yxy + B*x- 
where of particular interest is the coefficient 



D 



A X V — 



sin(4en)[[3 2 (t)-e 2 f3 2 x (t)] 



(18) 



which controls the entanglement between the 
two particles. As anticipated, it is equal to 
zero for n = and n = n cr . The other two A- 
coefficients are A xx = 



cos (2en) 
2^(t) 



M v = 



2 (2en) 



sin (2en) J_ 
2/31 (t) ^ 



sin 2 (2en) 2 j 
(t) e anCt 

We omit the explicit 



expressions for the coefficients B x and B^. 



At n = 



we obtain A%* — 



"2M 



= A^e 2 and 



Aq x ^z, which implies that, by the 
time when the collisions stop, the two wave packets ex- 
change the kinetic energies stored in the momentum dis- 
tributions in their respective rest frames. 

We finally mention, that the same results can also be 
obtained by representing both, the reflection from the 
wall and the two-particle collision as linear transforma- 
tions of the coefficients A xx , A™, A%", B x , B% and then 
multiplying the transformation matrices. 

In conclusion, we have described a new method for 
dealing analytically with interacting one-dimensional 
wave packets and applied it to the system depicted in 
Fig. [TJ The essence of the new method is a reduction of 
a quantum problem to classical calculations. The idea of 
this method should be generalizable to a broader class 
of problems with non-Gaussian shapes of the wave pack- 
ets, larger numbers of particles and different forms of the 
interaction potential. The unexpected result is that the 
entanglement between the two particles, that arises due 
to the collisions, vanishes by the time when the collisions 
stop. 
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SUPPLEMENTARY MATERIAL 

Note: In this supplementary material, we use the same 
notations as in the main article unless explicitly specified 
otherwise. 



cf>(— x;t), then consider the antisymmetrized superposi- 
tion of the original packet and the image, and, finally, 
take only the positive- a; part of this superposition: 



<f> a (x;t) = [(/>(x;t) 



-x;t)} Q(x). 



(S4) 



A. Full expression for the initial wave function in 
Eq.(2) of the main article 

The sign "~" in Eq.(2) of the main article means that 
the wave function needs to be slightly modified to exclude 
the mutual penetration of the two particles and the wall 
and then normalized. Since the tails are exponentially 
small, the details of this modification are not important. 
One just needs to introduce a cut-off factor, which varies 
sufficiently slowly as a function of x and y. A possible 
way of doing this is 



ip {c) (x,y) = Aexp^ 



[x - xmo} 2 



+ ip x ox - 



[y-yMo? 



x sin — 0(x)0(y — x) 

V v 



(SI) 



where N is a normalization constant, 0(x)0(y — x) are 
the Heaviside step functions. The sine- factor in Eq. lfSTj) 
preserves the continuity of the wave function at x = 
and x = y. 

Once the light particle approaches either the wall or 
the heavy particle, the behavior of the respective tail 
will evolve to follow the result of the "antisymmetrized 
image" procedure described in the next two sections. 



B. Reflection of a wave packet from the wall 

Here, we briefly summarize the "antisymmetrized im- 
age" method for calculating the reflection of a quantum 
wave packet from a hard wall. 

The Hamiltonian of a quantum system consisting of 
one particle and an impenetrable wall is 



H x — — 



if- 



2m x dx' 



+ S- G(-x) 



(S2) 



in the limit S —> oo. In the region x > 0, the parti- 
cle is free to move, but it cannot penetrate the region 
x < 0. Therefore, its wave function should continuously 
approach from the side of positive x and stay on the 
negative- a; side. 

Let us consider a wave function, which is well approx- 
imated far from the wall by a Gaussian wave packet 



(x;t) = N(t)exp 



IPxOX 



(S3) 



In order to respect the boundary condition at the wall, 
we introduce the mirror image of the above wave packet 



One can check directly that the above wave function is, 
in fact, an exact solution of the particle- wall problem. 

According to Eg. ((S3)) , the two wave packets <f>(x; t) and 
(f>{— x; t) have opposite group velocities. Initially, they 
move towards each other, then they meet and interfere 
at x — and, finally, separate again. Once they have 
separated, the main contribution to (jS4[) comes from 
the Gaussian (j>{— x;t). This is why the approximation 
(j> a (x;t) ~ — (f>{— x;t) is valid for sufficiently long times 
after the collision. 

To summarize, the wave function after the reflection 
from the wall is the wave function of the mirror image 
evolving in free space. What happens during the collision 
is easily computable, but we are not interested in it. 



C. Collision between two hard-core particles 

A collision between two hard-core particles is described 
by the Hamiltonian 



H 



+ 7p- +S-e(x-y) (S5) 



2m x 2m 



in the limit S — > oo. 

Before the collision, as long as the particles are suf- 
ficiently far from each other, the wave function of the 
system is well approximated as 



tp(x,y;t) ~ exp 



x M {t)Y 



x exp ■ 



2(3 x (t) 2 

[y-VM(t)Y 

2/3, (i) 2 



+ ipxox > (S6) 



We perform the standard transformation to the center- 
of-mass coordinate 



R 



m x x + m y y 



m x + m y 



and the relative coordinate 



r = x — y. 



(S7) 



(S8) 



The Hamiltonian in the new coordinates reads 
H = Hr + H r , where 



H R = 



2M'' 
Pi , 



(S9) 



H r = m + S ■ O(r) and S -> oo (S10) 
2/j, 



6 



with the total mass M and the reduced mass fi, respec- 
tively 



M 



m x + m y , 
m x m„ 



m,, 



(Sll) 
(S12) 



The dynamics of the center-of-mass coordinate R is just 
the free-particle evolution. At the same time, the Hamil- 
tonian H r is exactly the same as the one in Section \B\ 
The only difference from the one-particle situation is 
that now two variables are present. By analogy with 
the particle-wall collision, it can be shown that, in or- 
der to obtain the wave function after the collision, the 
initial wave function should be propagated in time as if 
the two particles were not interacting and then the sign 
of r should be changed but not the sign of R. This fact 
makes the result of this procedure non-trivial once we 
rewrite the wave function after the collision in the (x,y)- 
coordinates: 



J (m y - m x ) 2 2 

CX P { - n0 o,J, — ; — «ir 



(S13) 



2m 2 y 



x exp 



2p 2 (t)(m x +m y 

2m y x c (t) 
Pl(t)(m x +m l 



f3 2 {t)(m x +m y ) 2 



2m„ 



-y - Wxo 



m x + m y 



2m x (m v — m x ) 2m y (m v — m x ) 
x ex P i - o^,w \o x y + o->rZ, v> x V 



Pl(t)(m x +m y ) 2 Pl(t)(m x + m y ) 



x exp 



x exp 



2m? 



P 2 (t)(m x +m y 



(m y -m x ) 2 ^ 2 



2(3 2 (t)(m x + m y ) 



(m y - m x )x c (t) 



x + ip x0 - 



by 1 1 L X 



xl(t) 



Pl(t)(m x + m y ) m x +m y 2(3% J' 

where x 2 is given in the main article. The momenta after 
the collision PxO m l+ v my and Pxa Zl+Zl appearing in 
Eq. (|S13I) follow from the classical energy and momentum 
conservation. 



D. Probability distributions for the classical 
multiple collisions problem 

As a part of the full quantum solution, we need to 
solve a classical problem with the following initial con- 
ditions: A wall is at the origin. The light particle has 
the initial momentum p x o and a definite initial position 
xmo, i«e. the initial probability distribution of positions 
is Px,o(x m ;0) — S(x m — xmo)- The heavy particle has 
zero initial momentum and the Gaussian distribution of 
possible initial positions 



P y ,o(y m ;0) = exp < - 



[Vm - 3/mo] 



2Ac 



(S14) 



y,0 



Below we calculate the evolutions of the marginal prob- 
ability distributions P K . n (x m ;t) and P y ,n(ym',t) as func- 
tions of the number of collisions n occurred. In order to 



do this, we first fix an arbitrary initial position of the 
heavy particle y m o and solve the problem in this specific 
case. Then we address the evolution of the probability 
distributions. 



1. Solution for a given y m o 

At t = 0, we have the following initial situation: the 
light particle is at xmo with a momentum p X Q directed 
towards the heavy particle, which is at the position y m Q. 

Before starting with the quantitative analysis, let us 
summarize the qualitative expectations. Due to the small 
but nevertheless non-zero mass ratio e = y/m x /m y <C 1, 
the light particle will lose its momentum and accelerate 
the heavy particle at each collision. This transfer of mo- 
mentum lasts as long as the light particle is faster than 
the heavy one. After the maximal number of collisions 
"max, the light particle will become too slow to catch up 
with the heavy one. Because of the smallness of e and the 
conservation of energy, the velocity of the heavy particle 
after the last collision is much smaller than the initial 
velocity of the light particle v x q = Pxo/m-x- 

To compute the position of the heavy 
particle, we introduce the function 

n(y m a',t) = n, which represents the number of col- 
lisions before time t. With such a definition, 



2/m(2/mo; t) = y(n) + [t- i(n)] v y (n), 



(S15) 



where y(n) and t(n) denote, respectively, the position 
and the time of the n-th collision and v y (n) the velocity 
of the heavy particle after the n-th collision. 

The velocities v y (n) and v x (n) after n collisions can be 
obtained using the following iteration relations 



v x (n + 1) = 
v y (n + l) = 



x i \ 
—v x (n) 



2m 



V—v y (nUS16) 



2m x m y -m x W01 ^ 
-v x (n) + Vy(n),(S17) 



in. 



which lead to 



v x (n) = v x0 cos(nip), 
v y (n) = v x0 esin(n(f), 



where 



ip = arctan 



2c 



2e. 



(S18) 
(S19) 



(S20) 



The maximal number of collisions can be determined 
from the condition 



which gives 



arctan (-) ir 
~ 4e' 



arctan 



(A) 



(S21) 



(S22) 
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To calculate the positions and the times of the colli- 
sions, another set of iteration relations can be obtained: 



v x (n)+v y (n) 
y(n + 1) = — r-r r^y{n), 



t(n + 1) - i(n) 



2 



y(n). 



(S23) 
(S24) 



Expanding expression (jS23[) in powers of the small pa- 
rameter e up to the second order, we obtain 



y(n + 1) « (l + 4ne 2 ) y(n). 



(S25) 



Since we expect the motion of the heavy particle to be 
slow, the following approximation is valid 



y(n + 1) - y(n) dy(n) 



1 dn 
and the solution to this equation is 

y(n) = 2/moexp(2n 2 e 2 ). 



Ane 2 y{n) 



(S26) 



(S27) 



Using Eq. (IS27|) and making the same approximation, we 
obtain from Eq.fS24|) 



t(n) 



'Ml 



ml) 



V X 



1 - e 2 -n 2 



n+ 3 



(S28) 



Important to know is the inverse of this function, which 
gives us the number of collisions occurred during time t 



n(y m o;t) 



tv x0 tv xQ {2tv x0 + y m0 ) (tv x0 - 2y m0 ) 



2y 



mO 



I2y 



mO 



(S29) 

The right-hand-side of the above equation is implied to 
be rounded from below. The first term in Eq. (|S29l) is 
what we expect from an infinite mass particle located at 
y m Q. The minus sign in front of the second term is also 
expected because a non-zero mass ratio e means that the 
same number of collisions takes longer time. 

The position of the heavy particle between the n-th 
and the n + 1-th collisions is 



Vm{ymo;t) = y(n) + [t- i(n)] v y (n). 



(S30) 



Likewise, the position of the light particle x m (y m o; i) be- 
tween the two collisions can be written as 



x m (ymo;t) = \y(n) - [t- t{n)] v x (n)\ 



(S31) 



2. Evolution of the classical probability distributions 

After obtaining the results for a given initial position 
of the heavy particle, the evolutions of the marginal pro- 
bability distributions P x ^ n (x m ;t) and P y , n {ym',t) can be 
calculated. These probability distributions refer to the 
moments of time between the collisions, when, given the 
assumption of the narrow wave packets, the number of 
the preceding collisions can be assumed to be the same 



for any value of y m o. We now focus on the shapes and the 
widths of these distributions as functions of the number 
of collisions n. 

The evolution of the probability distributions can be 
divided into two alternating stages: the evolution be- 
tween the collisions and the change due to the collisions. 
During the first stage, i.e. for a fixed value of n, the 
probability distributions P Xl n(x m ',t) and P y ,n(ym',t) only 
shift along the x m - and y m -axes with respective velocities 
Pxn/mx and p yn /m y . Therefore, the shapes and widths 
of the distributions only change during the second stage 
when n changes. 

As mentioned in the main article, two trends in the evo- 



lution of P x . 



;t) and -P yi n(y m ;i) can be anticipated 



qualitatively, namely, an initial broadening of the dis- 
tribution P x , n {xm',t) for the light particle and an initial 
narrowing of the distribution Py t n(y-m]t) for the heavy 
particle. Both trends can be understood with the help 
of Fig. IS 1 1 which represents the time evolutions starting 
from two different initial conditions: Case A and Case B. 

In the beginning (i = ii), the light particle is at the 
same position in the both Cases but the heavy particle 
is not. This is why at a later time (t = t2), the first 
collision has already happened in Case A but not yet in 
Case B. In Case B, the particles finally collide at t = 
t^. After one collision in the both Cases (t = £4), the 
two velocities of the light particles are the same and the 
two velocities of the heavy particles are also the same. 
However, during the time delay between t-2. and £3 the 
heavy particle was moving faster in Case A. This is why 
the difference in the heavy particle's positions in the two 
Cases decreases (distance between the red dashed lines). 
On the contrary, the distance between the two positions 
of the light particle has increased during the same interval 
from zero to a finite value (distance between the green 
dashed lines). 

The changes of the relative distances in Fig. [STI 
are exaggerated. In reality, these changes are rather 
small. However, they accumulate after multiple colli- 
sions. Moreover, as shown below, the two trends de- 
scribed above additionally enhance each other. 

We now calculate the narrowing of the distribution 
Py,n(ym',t). For this purpose, we investigate the evo- 
lution of the distance between two different initial posi- 
tions of the heavy particle. Denoting the initial distance 
as (Ay)o and the distance after n collisions as (Ay)„ we 
obtain 



(Ay)„_i - (Aj/)„ = Av n At n , 



(S32) 



where Av n is the difference in velocity of the heavy par- 
ticles during the time delay. It follows from Eq. (|S19p 
that 



Av n = [sfn(n<£) - sin((n - l)tp)]ev x0 . 
The time delay in Eq. ([S32)l is 

(Ax)„_i + (Ay)„_i 



AU, 



[cos((n -!)</?)-£ sin((n - l)^)]^ ' 



(S33) 



(S34) 
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t = U 



t = u 



t = ti 



Case A: '4 



Case B: 



_QQ_ 



XXX 



i 



FIG. SI. In this sketch, the two rows represent two different initial conditions. The upper and the lower lines of figures 
represent respectively Cases A and B described in the text. The columns represent the evolution in time. The arrows depict 
the velocities of the particles. 



In the numerator, (Ax) n _i is the distance between the 
light particles after n — 1 collisions. The denominator of 
Eq. (|S34l) is the difference between the velocities of the 
heavy and the light particles. 

The distance between the light particles can be ob- 
tained from 

(Ax) n = (Ay) n _i + At n [cos(n^) + esin((n - l)ip)] v x0 . 

(S35) 

In the context of Fig. IS1[ the second term in Eq.fS35j 
corresponds to the distance the light particle travels in 
Case A between t = t% and t = t%. Inserting (|S34j) into 
(|S35|) we obtain 

n-i 

(H « (Aaj)„_i + 2(Ay) n _i « £ 2(Ay) i . (S36) 



Substituting (fS33|) and (|S34|) into (|S32]) and expanding 
the result in terms of e up to the second order, we further 
obtain for a large n 



(Ay) n+1 - (Ay)„ 



-4e 2 £(Aj/), 

2 = 



(S37) 



which, in the continuous limit, can be rewritten as 
d(Ay) n 



dn 



-4e' I {Ay) n , dn' . 



(S38) 



The solution of this equation is 

(Ay)„ = (Ay) cos(2en). 



(S39) 



The consequence of the dependence (|S39[) is that the 
shape of P y ^ n {y m \t) remains Gaussian. This originates 
in the fact that the scaling factor cos(2en) in Eq. (|S39[) 
does not depend on y m o, i.e. the distance between two 
arbitrary points of the distribution P yi n(ym]t) becomes 
rescaled by the same factor as long as all points of the 
distribution have experienced the same number of colli- 
sions. Therefore, it follows from Eq. (|S39|) that the width 
of Py,n{ym', t) evolves as 



(Aa y ) n = A<7 y \ cos(2en)|. 



(S40) 



The explicit form of P y ,n(2/m; t) is thus 



Py,n(y m ]t) = CXp 



[Vm - S/mCQ] 

(Aa y )n 



(S41) 



where y M (t) = 2/ m (j/Mo; t) is given in E q.(|S30|) . 

To calculate (Ax) n , we substitute (|S39|) into (|S36I) . 
which gives 



(Ax) n «2(At/) ^cos(2ei), 



(S42) 



i=0 



and then perform the summation with the help of the 
following approximation: 



n-l 



cos(2ei) 



i=0 



1 



cos(2ea;)c?a; = — sin(2e(n — 1)). 



(S43) 

Substituting Eq. (|S43p into (|S42I) and approximating 
n-l«B,we obtain 



{Ax) n w ^^sin(2en). 



(S44) 



Since the scaling factor sln (^ e ") j n Eq. (|S44[) does not 
depend on y m Q and Py in (y m ;i) always has the Gaus- 
sian shape, the distribution P x , n {x m ',t) also acquires the 
Gaussian shape. Its width (Aa x ) n evolves as 



{A<T X 



| sin(2en)|. 



(S45) 



Equations (12) and (13) of the main article for the evo- 
lution of the particles' coordinates follow from Eqs. (|S39p 
and (1S44[) and from the definitions Xyi(t) = x m {yyio;t) 
[the RHS given by Eq.([S31J] and y M {t) = y m (yuo\t) [the 
RHS given by Ea. (jS30)) ]. 



E. Ballistic spread of the wave packets and the 
validity of the narrow-wave-packet assumption 

The validity of our approximation is limited by the 
narrow- wave-packet assumption that both |/3a;(t)| and 
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should be much smaller than yyi(t), which for the 
sake of an estimate can be further approximated as umo- 
We limit ourselves to the case |Ac(i)| 3> \f3 y (t)\. There- 
fore, the critical condition invalidating our approxima- 
tion is 



Vmo = \P*(t)\ 2 = ° 2 0x 1 + 



t 2 h 2 

1 4 



t 2 h 2 

11 2 f-r '2 



(S46) 



Denoting the number of collisions required to reach this 
condition as ribai, then approximating the corresponding 
time as t « Sj/Mo^bai/^o and, finally, substituting it in 
Eq. (|S46|l . we obtain 



"-bal 



m x a 0x v x o 
Ah 



(S47) 



The condition for the validity of our approximation be- 



< "bal, 



(S48) 



where, as derived in the previous section, 
Condition (|S48[) is thus equivalent to 



e m x a 0x v x0 

Tit) 



> 1. 



(S49) 



Let us now give a concrete example when the left-hand- 
side of inequality (|S49|) is approximately equal to 1 . This 
critical situation may be realized when the light particle 
is a big molecule with m x = 10 _20 g and the heavy parti- 
cle is a very small dust grain of mass m y = 10 _10 g. The 
initial velocity of the molecule is v x o = 10 4 ™. The initial 
wave packet widths of the two particles are ao x — 10~ 6 cm 
and a 0y = 10~ 10 cm. 



In order to satisfy condition (|S49|) starting from the 
above numbers, one should increase the value of either of 
the parameters appearing in the product e rn x a 0x v x o. 



